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by describing the efficient organization of amplitudes resulting from the conjectured 
duality between color and kinematic structures on constituent graphs. 
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1. Introduction 

Over the past few years there has been spectacular progress [H,[2}[3jlll|5l|6l[71IEJ 
[H [101 HH [12] in writing down the integrands of loop-level scattering amplitudes of 
the M = 4 super Yang-Mills (sYM) theory [13] in the planar leading-color sector, 
or 't Hooft limit. Similarly, for the integrated planar multiloop amplitudes there 
are recent important advances, notably the all-loop resummation of four and five- 
particle amplitudes proposed by Anastasiou, Bern, Dixon, Kosower, and Smirnov [14J, 
and by the fully analytic evaluation of the six-particle two-loop maximally-helicity- 
violating amplitude pj)]; beyond that, there are many hints, both at weak and strong 
coupling [161 HZl UHl H9l EQl ETJ [22l [231 El], that the planar sector may soon be solved. 

The subleading color contributions to N = 4 sYM amplitudes are less well studied— 
save for the one loop case, where all subleading-color amplitudes may be obtained from 
the planar sector ones (see e.g. [H [25]). Beyond one loop, the literature contains only 
the full non-planar contributions of the four-point amplitude, worked out at two [26] . 
three [271 123 [29] and four [30] loops at the level of the integrand. Notably, these four- 
point amplitudes also happen to be some of the few amplitudes known to be valid 
in higher dimensions D > 4 [301 EI] in the dimensionally oxidized M = 4 theory 
(dimensional reductions of the D = 10 sYM theory). 

In this contribution to the review, we will discuss some generally applicable methods 
for constructing multiloop multiparticle integrands of amplitudes, valid for TV = 4 sYM 
in four and higher dimensions, including non-planar contributions. Much of what is 
covered also applies to less supersymmetric theories, including pure QCD; however, the 
focus and examples in this chapter will be on the M = 4 sYM theory. We will exhibit 
the general organization of unitarity cuts and amplitudes, specifically the maximal 
cut method [H [281 EQ] and the imposition of a duality between color and kinematic 
structures of amplitudes [321 [29]. Interestingly, the duality allows one to recognize 
that all contributing terms of the amplitude are fully specified by some small set of 
independent "master graphs". A discussion which only applies to M = 4 sYM is the 
review of special recursive cuts known as "two-particle cuts" [261 E31 [30], and "box 
cuts" [30], relying on the fact that all four-point amplitudes in M = 4 sYM factorize 
into a product of the corresponding tree amplitude times a universal function. 

Loop amplitudes of massless theories in four dimensions suffer from infrared 
divergences (and frequently also ultraviolet divergences). A standard approach for 
avoiding these infrared divergences is to dimensionally regulate the amplitudes, that 
is, to work in D = 4 — 2e with e < 0. In this scheme, it is thus critical to compute loop 
amplitudes in D > 4 dimensions. The methods presented in this review are well suited 
for this task. 

Throughout this contribution we make use of the unitarity method [M] and 
generalized unitarity [351 ESI E7J , which we consequently will briefly outline. However, 
deeper understanding of details, as well as applications of this approach, may be found 
in some of the parallel chapters of this review, e.g. Bern and Huang [38], Britto [39] , 
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and Ita [40] . 

The organization is as follows. In section [21 we set up notation and organization of 
multiloop integrands and cuts. In section [3], the method of maximal cuts is discussed in 
some detail, focusing on both the systematic identification of needed cuts as well as the 
computation of these using tree amplitudes as input. Complimentary to the method 
of maximal cuts, in section HI we review two particularly simple classes of multiloop 
cuts, relevant to the M = 4 sYM theory: two-particle cuts, and box cuts. Finally, in 
section [51 we discuss the set of powerful constraints that one can impose on multiloop 
gauge theory amplitudes by considering the color-kinematics duality. 

2. Preliminaries 

2.1. Integral representations for amplitudes 

To systematically employ the unitarity method [34" | |35 | [36], it is convenient to organize 
the amplitude using a sufficiently general integral representation. We note that one 
can express any local massless gauge theory amplitude in D dimensions in terms of 
local interactions and non-local propagators ~ 1/p 2 , regardless of the particle content 
(cf. Feynman gauge propagators and vertices). With this in mind, one can write any 
amplitude as 

d LD i local 

where d LD i = Yl^j=i d D @j is the usual Feynman diagram integral measure of L 
independent D-dimensional loop momenta £j, and T is the set of all L-loop graphs, 
with n external legs, counting all relabeling of external legs. Corresponding to each 
internal line (edge) of the i th graph, we associate a propagator 1 / p? , which is a function 
of the independent internal and external momenta, £j and kj, respectively. The local 
numerator functions, here only schematically indicated, include information about the 
color, kinematics and states, etc. This representation thus mimics the behavior of 
Feynman diagrams, but as we will see the local numerators need not be calculated 
using off-shell Feynman rules. Instead it is often more efficient to indirectly infer these 
numerator functions by imposing general constraints and by matching to physical on- 
shell information gleaned from generalized unitarity cuts. 

Clearly the way of expressing an amplitude in (JT)) is not unique. One can always add 
more propagators to any graph by multiplying by unity, expressed as pi/ pi, absorbing 
the numerator p^ into the graph's local numerator factor without changing its locality. 
To remove this freedom, we can insist that the graphs we use have only cubic (trivalent) 
interactions. That is, we can hide the contact terms that one usually encounters in 
gauge theories through the repeated multiplication of pl/pt factors until every graph 
has the same number of propagators as a cubic graph: n + 3L — 3. Thus, we use a more 
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constrained representation of the amplitude, 

71 ~ 2-^ / (Ott\LD C. „2 2 2 . 2 ' I 2 ) 
i6 r 3 ^ ^ Z?r ' ^ PiiPi2Pi 3 Pim 

where T 3 is the set of all purely cubic graphs, and m = n + 3L — 3 for all graphs. 
Furthermore, we have, for convenience, separated out the kinematic parts Ni and the 
color parts Cj of the local numerator factor, and normalized the numerators by dividing 
out by the standard (Bose) symmetry factor Si of each graph. 

One might well wonder if it is always possible to factorize the kinematic and color 
terms for each graph. Indeed, this requirement imposes severe constraints on the local 
numerator. Nonetheless, when all particles are in the adjoint representation of the 
gauge group, as is the situation for the theories we will consider, this factorization is 
possible [32] . It can be seen from the Feynman rules of gauge theory: if all particles are in 
the adjoint the only color structure that can appear in cubic vertices are the gauge group 
structure constants f abc . Hence, every purely cubic Feynman diagram has a unique 
color factor associated to it. For the quartic contact vertices of gauge theory, there are 
three different color structures: f aiCL2b f ba ^ a * and its permutations. If we promote such 
a contact term to be a part of a cubic graph using f aiCL2b f ba s a i _s. p 2 2 /pl 2 f aia2b f ba3a4 , 
then this term is proportional to the same color factor as the non-contact contribution 
in the p\ 2 = (pi + P2) 2 channel. Thus, by this construction, the color factor of each 
cubic graph in (j2J) is unique, up to a normalization factor. They are given by 

d = Yl j u , (3) 

j€Vi 

where the product runs over all the cubic vertices in the i th graph. Here, we use structure 
constants f abc with a normalization differing from normal conventions, f abc = i\^2f abc = 
Tr([T a , T b ]T c ) , and with hermitian generators T a normalized via Tr(T a T fe ) = 5 ab . 

The kinematic local numerators JVj are by definition polynomials in the loop 
momenta but may be taken to be either polynomials or rational functions of external 
kinematic variables depending on the convenience of the situation at hand. Needless to 
say, Ni also depends on the polarizations, helicities and spinors, etc, according to the 
details of the external states. The maximal degree of the polynomial in loop momenta 
can be determined by dimensional analysis: an n-particle gauge theory amplitude in 
D = 4 has dimension 4 — n; thus, after accounting for the dimension of the integral 
measure and propagators, Ni has dimension n + 2L — 2. This dimensionality coincides 
with the number of cubic vertices in each graph, as is expected of a gauge theory where 
each cubic vertex has one power of momenta. The dimension of iVj then has to be 
distributed between the factors of external and internal momenta. This means that in 
a generic gauge theory, one would expect that the loop momentum dependence in Ni 
comes in monomials of the form 

lj\f,...,£ n+2L - 3 ,£ n+2L ~ 2 G Ni, (4) 

where i v is a mnemonic for any Lorentz tensor constructed out of p factors of loop 
momenta i^. For special gauge theories, the highest degree monomials are absent for 
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all iVj. This property can be attributed to the good ultraviolet behavior of the theory, 
as the ultraviolet regime is indeed the corner of phase space in which £^—7-00. 

Although the ultraviolet behavior of a theory is interesting in itself, for the purposes 
of this review we will only note that such behavior constrains the form of the amplitude. 
In the absence of other constraints, the fewer powers of £ that one must consider, the 
smaller the space of possible numerator polynomials will be. This can greatly facilitate 
amplitude construction. The set of well-known gauge theories with good ultraviolet 
behavior includes many supersymmetric Yang-Mills theories. In particular, maximally 
supersymmetric M = 4 sYM is known to reduce the maximum power of the monomials 
(jlj) by at least four relative to pure Yang- Mills theory jJH 02} [331 S3]- Theories with 
A/"-fold supersymmetry are expected to reduce the maximal power of i by at least M 
in ([4]), see e.g. [44] . In some cases, the cancellations of the high powers of i are even 
stronger than this naive counting (321 [301 HH] • 

For completeness, we also mention an alternative integral representation frequently 
used in the literature. For amplitudes with many external legs (compared to the number 
of spacetime dimensions), it can be cumbersome to use an integral representation with 
only cubic vertices. For example, the standard one-loop integral basis in four dimensions 
consists of box, triangle and bubble integrals; no higher polygon integrals are needed 
for massless D = 4 theories (see e.g. [39J). A similar basis is expected to exist at the 
two-loop level jlEj . To allow for such cases, we use a sufficiently general (over-complete) 
integral representation 

A n L) = Yl Gk Ri ' k [ (2n) LD v 2 v 2 v 2 ■ --v 2 ' ^ 

where the first sum runs over color factors Gk in the trace basis (TB), of say SU(N C ), and 
the second sum runs over n-point L-loop graphs as in (CQ). The Ri t k are rational functions 
of external momenta, and are polynomials of the independent loop momenta. As 
desired, one can generally restrict the number of propagators m < \D~\L ([•] is the 
ceiling function), because integrals with m > \D~\L are reducible {e.g. using Melrose, 
Van-Neerven-Vermasseren, or Bern-Dixon-Kosower reduction and similar methods |47j ) . 
although the exact upper bound depends critically on the subtleties of dimensional 
regularization at each loop order. Similar to the representation (T5]) the ultraviolet 
behavior of the theory under consideration imposes useful constraints on the form of 
Pi,k- For example, the good ultraviolet behavior of M = 4 sYM can be used to show 
that only scalar boxes contribute at one loop in D = 4 pQ; thus, one can set P^ = 1 for 
all boxes and Pj^ = for all other one-loop integrals in this theory. 

As a recent alternative to the mentioned general integral representations, we remark 
on an interesting new set of chiral integrals that have been introduced specifically 
for planar M = 4 sYM using momentum twistor space notation [TT] . Indeed, using 
special properties of a theory, such as dual (super-) conformal symmetry [161 EH] 
(or Yangian structure [20]) of M = 4 sYM, one may construct more constrained 
integral representations, which consequently increases the power of any amplitude 
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calculation [3JI11EJE]- m similar spirit, the duality between color and kinematics [3"2|I2§] 
discussed in section [51 is conjectured to provide an constraining organizing principle for 
the integrals in more general gauge theories. 

We note that the color-kinematics duality will critically rely on the cubic integral 
representation (J2J), whereas in sections [3] and HI which treat generalized unitarity cuts, 
either integral representation ([2]) or ([5]) is applicable. 

2.2. Generalized unitarity cuts 

Unitarity of the S"-matrix requires that loop amplitudes must satisfy certain properties 
relating them to sums over products of lower-loop amplitudes. The goal of the unitarity 
method [Ml 1351 136] is to construct a compact expression that obeys all constraints 
required by unitarity for a particular amplitude. For massless theories it is sufficient to 
satisfy all generalized unitarity cuts in D dimensions, this guarantees that no logarithmic 
or rational terms can go undetected. Strictly speaking only a subset of these, a " spanning 
set" of cuts [30J, are necessary to guarantee the satisfaction of all other cuts. 

To identify contributions to the amplitude, we will compute generalized unitarity 
cuts on the level of the integrand, expressing them as a product of on-shell 
subamplitudes, 



where Au\ are the subamplitudes, which may be either tree or lower-loop amplitudes. 
The number of cut lines for each cut (the number of on-shell pj = internal legs) is 
denoted by c. 

The cuts can be computed directly from the theory by feeding in the corresponding 
subamplitudes and summing over the intermediate states, as in (Q. They may, of 
course, also be computed from graph-based integrand representations of the full loop 
amplitude, (T5]) or ([5]), when available. In this latter case, the sum over intermediate 
states has already been carried out, and the generalized cut can be formally obtained 
by replacing the propagators of the cut lines by delta functions, 



and collecting the terms with the most delta functions. Since the full integrand of the 
amplitude under consideration is typically a priori unknown, the main purpose of these 
types of cuts is to apply them to an ansatz of the amplitude. 

The unitarity method instructs us to reconstruct amplitudes using the information 
from the set of all generalized unitarity cuts of a theory. This set is overcomplete; 
consequently, there are different strategic approaches on how to most efficiently extract 
the relevant information for building the amplitudes. For example, one might well 
consider a strategy of using the fewest number of cuts, with each individual cut thereby 
providing the most information about a given amplitude. The most extreme such 
situation is exemplified by the one-particle cut of an n-point L-loop amplitude, relating 




(6) 



states 




(7) 
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it to a (n + 2)-point (L — l)-loop amplitude sewn with itself. This cut has the minimal 
number of cut propagators, and in principle contains the full information of all the other 
generalized cuts of the L-loop amplitude (ignoring any subtleties with possible collinear 
and soft divergences of the cut line). Without an organizing principle, disentangling 
this information can be prohibitive. In the planar-limit of M = 4 sYM, dual conformal 
symmetry (or Yangian invariance) appears to provide such an organizing principle [11] . 
See also [48] for a general discussion of single-line cuts at the one-loop level. 

In the next section, we consider a different strategy, using the method of maximal 
cuts. Rather than maximizing the information coming from each cut, we maximize 
the number of cuts and similarly the number of cut lines, ensuring that each cut 
provides as small an amount of information as possible, thus facilitating the piece- 
by-piece construction of the amplitude. 

3. Maximal cuts in D = 4 and general dimension 

The method of maximal cuts [U EHl [30] offers a particularly efficient means for 
determining multiloop amplitudes. In this method, we start from those generalized 
cuts that have the maximum number of cut propagators (maximal cuts). These cuts 
provide an initial ansatz for the amplitude. Next we systematically correct the ansatz 
by considering the information provided by near-maximal cuts. These sequentially 
reduce the number of cut propagators one by one. In this process, all potential contact 
contributions are determined. 

Generalized unitarity has a wide range of applicability: amplitudes of generic 
theories in generic spacetime dimension can be constructed. Since maximal and near- 
maximal cuts are special cases of generalized unitarity cuts, they can be used in 
many theories and space-time dimensions. However, there are some mild assumptions 
necessary if one wishes to organize around cubic vertices and cut all propagators: 
it is important that massless on-shell three-point amplitudes are non-vanishing and 
non-singular [37J |49], for appropriate choices of complex cut loop momenta [501 El]- 
Moreover, it is important that the space-time is of dimension D > 4, in order to allow 
for one to impose both momentum conservation and on-shell conditions of the legs 
meeting at a three-point vertex without having collinear momenta^. 

The strategy of maximizing the number of cut propagators is known to be a powerful 
method for computing one-loop M = 4 sYM amplitudes in four dimensions, because 
only box integrals appear pQ. As observed by Britto, Cachazo and Feng [37], taking a 
quadruple cut, where all four propagators in a box integral are cut, freezes the four- 
dimensional loop integration. This allows the kinematic coefficient of the box to be 
straightforwardly computed directly in terms of the cut. The use of complex momenta, 
as suggested by twistor space theories [51], makes it possible to define massless three- 

| These assumptions are needed if one wants to cut down to cubic vertices. In theories that have no 
on-shell three point vertices amplitudes may instead be organized around quartic vertices, which has 
no such restrictions [52] . 
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vertices and thereby to use quadruple cuts to determine the coefficients of all box 
integrals including those with massless external legs. 

The idea behind the quadruple cuts has been generalized by Buchbinder and 
Cachazo [49] to two loops using hepta- and octa-cuts of the double box integral 
topology. Although the two-loop four-point double-box integral only has seven physical 
propagators, it secretly enforces an additional spurious eighth constraint, yielding an 
octa-cut which localizes the integration completely. These types of cuts are now known 
as a part of the "leading-singularity" technique, which is similar to the maximal 
cut method in that it is based on cutting a maximal or near-maximal number of 
propagators (53J Ej, but in addition it may make use of extra conditions from spurious 
singularities that are special to four dimensions. 

The strength of the method of maximal cuts has most clearly been seen in 
amplitudes that have a large number of multiloop integrals. For example, the four- 
point amplitudes of M = 4 sYM have been evaluated using maximal and near-maximal 
cuts at five loops [4] in the planar sector, and similarly for the full non-planar amplitudes 
at three [28] and four loops [30], the latter containing 50 distinct integrals, each with 
multiple terms. 

3.1. Working definition of maximal and near-maximal cuts 

To make the concept of maximal cuts clear, we will introduce some proper definitions. 
First we note that one can represent generalized cuts ([6]) as graphs: vertices (nodes) 
representing on-shell subamplitudes, and lines (edges) representing on-shell propagators. 
Maximal cuts are simply those generalized cuts (graphs) that have the most possible on- 
shell cut conditions, given three numbers: n the number of external legs, L the number 
of loops and D the number of dimensions of the loop momenta. For any generalized 
cut, there are two competing upper limits to the number of cut conditions c that one 
can impose: 

c < n + 3L — 3 , (cubic graph limit) 

c < DL , (loop momentum freeze-out) (8) 

where the first limit comes from counting the internal lines of a cubic n-point L-loop 
graph, and the second by counting the degrees of freedom of the independent loop 
momenta. It is clear that a cut consisting of only cubic vertices can have no more 
physical singularities as the integrand is completely local. Similarly, a cut where all the 
degrees of freedom of the loop momenta are frozen can have no more singularities in the 
internal phase space. 

We define maximal cuts to be those generalized cuts which saturate one of the 
above bounds^], that is, 

cmc — min{n + 3L — 3, DL} . (9) 

§ In some unitarity cuts, the "loop momentum freeze-out" bound, which counts the total degrees of 
freedom of the loop momenta, is not a good measure of what happens locally in each subgraph. For 
example, a direct product of a one-loop triple cut and a one-loop penta-cut would be a two-loop maximal 
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For example, for D = 4 and L = 1 the maximal cut is equal to the quadruple cut 
cmc — 4. For n = 4 and L = 2 it is the hepta-cut cmc = 7, etc. (For comparison, we 
note that if we work in four dimensions, then for loop orders L < n — 3 the maximal 
cuts coincide with the leading singularity cuts as both have 4L cut lines.) 

To define the set of near-maximal cuts, we consider those cuts with fewer cut lines 
than the maximal cuts c < Cmc- m particular, a next-to-maximal cut has c = cmc~ 1 cu t 
lines, and a next-to-next-to- maximal cut have c = cmc ~~ 2 cuts, and so on. For example, 
in this classification, we refer to both one-loop triple cuts in D = 4 and quadruple cuts 
in D = 5 as next-to-maximal cuts (for n > 4). 

3.2. An algorithm for the maximal cut method 

We give here a general outline for constructing the loop amplitudes using the method 
of maximal cuts. In words, the algorithm is as follows: 

(i) Enumerate all possible maximal cuts given the multiplicity n, loop order L, and 
dimension D. For example, if n < D, or more generally L > in — 3)/(-D — 3), the 
maximal cuts are in one-to-one correspondence with the cubic graphs of loop order 
L and multiplicity n (as follows from iQ). 

(ii) Compute the maximal cuts at the level of the integrand (i.e. ([6])). If the maximal 
cuts do not freeze the loop momenta, make sure to capture this residual dependence 
using fully Lorentz covariant variables of the loop momenta (Lorentz products or 
Levi-Civita tensor contractions). 

(iii) Using the information from the maximal cuts, construct an initial ansatz for the 
amplitude. By the Lorentz covariant construction, the expressions obtained from 
the maximal cuts are automatically pushed to off-shell loop momenta (the validity 
of any such continuation will be ensured by the consecutive near-maximal cuts). 

(iv) Enumerate the next-to-maximal cuts. Calculate each cut using ([6]), then compare 
to the result of each cut as obtained from the amplitude ansatz constructed in the 
previous step. If the difference between the true cut and the cut of the ansatz is 
nonzero, then correct the ansatz by incorporating said difference. 

(v) Continue iterating the previous step for the next-to-next-to-maximal cuts, and so 
on; each step with one fewer cut condition imposed. 

The iteration can be stopped at any point when the amplitude ansatz appears 
to receive no more corrections. A proof that the ansatz is correct will follow from 
evaluating it on a complete set of "spanning cuts" [30], which are cuts that are known 
to fully determine n-point L-loop amplitudes of any theory (see [38] for more details). 
Alternatively, the iteration may be continued by removing cut conditions one by one 

cut in four dimensions, since it is an octa-cut. However, since the one-loop penta-cut does not exist in 
four dimensions, one could argue that it should not be counted. Nonetheless, in order to simplify our 
definition ([9]), we simply count such situations as valid cuts that happens to evaluate to zero. 
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until the iteration naturally terminates by generating its own set of complete "spanning 
cuts" . 

If we represent the amplitude using the integrals in (j2J), then the amplitude is 
entirely specified by the integral numerators iVj. The idea of the above algorithm is 
then that the numerators can be expanded ajJJ 

N l = Nf AC + N? MC + Nf MC + ... , (10) 

where Nf 10 is the initial ansatz determined by the maximal cuts, and the corrections 
from the near-maximal cuts take the schematic form 

m m 
j=l j,k=l 

where the p 2 . are inverse propagators and P^j and Pi-j t k are polynomials fixed by the 
next-to- and next-to-next-to-maximal cuts, respectively. The contributions from next- 
to-maximal cuts must be multiplied by at least one inverse propagator since these 
contributions must vanish on maximal cuts where all pi. = 0. Similarly, the corrections 
from the next-to-next-to-maximal cuts be multiplied by at least two inverse propagators 
since they vanish on all next-to-maximal cuts, and so on. Furthermore, since the inverse 
propagators carry dimension, the dimension (or degree) of the local polynomials must 
reduce for each iteration step, yet by locality remain non-negative. This is why the 
algorithm naturally terminates. 

In the following sections, we will review how to compute maximal cuts using gauge 
theory tree amplitudes as input. 



3.3. Cuts with on-shell three-vertices 

In this part, we direct our attention to the evaluation of the maximal cuts. First, we 
will consider maximal cuts consisting of only three-point tree amplitudes 

•c \ 1 a tree a tree a tree a tree /10\ 

1 3 ,(1) 3 .( 2 ) 3 .(3) ' ' ' 3,(m) > \ lZ ) 

states 

where the number of cut conditions is given by c = n + 2L — 2. Many of the interesting 
features of this cut are exposed by studying the kinematics of a single on-shell tree A^ ee . 

Consider three massless momenta p, q and r = — p — q meeting at a three-point 
vertex, as in figure [H one has the following constraints: 

p 2 = q 2 = p ■ q = , (13) 

where the orthogonality p ■ q = follows from the masslessness of the third leg 
r 2 = (p + q) 2 = 2p ■ q = 0. For real momenta in Minkowski signature (1,D — 1), these 

|| It should be noted that this expansion is not unique. For example, contributions of the amplitude 
corresponding to contact terms can be shuffled between different graphs without affecting the amplitude, 
nor the generalized unitarity cuts. The algorithm of the maximal cut method simply gives a 
representation of the amplitude that satisfies all unitarity cuts; thus, much of the reparametrization 
freedom of the full amplitude is present in this representation. 
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equations have only degenerate solutions where the momenta are collinear, p || q. For 
nondegenerate solutions we may either work in (2, D—2) signature, or more generally, we 
may explicitly complexify the momentum space (the (2, D — 2) signature can be thought 
of as a special case of complexfied momenta). For example, an explicit nondegenerate 
solution can, choosing a convenient Lorentz frame, take the tag 

p»= (1,1,0,0,0,. ..,0), g" = (0,0,l,i,0,...,0). (14) 

Interestingly the solutions in D < 4 are strictly degenerate; thus the on-shell three- vertex 
is only supported by kinematics in D > 4. The case D = 4 is the most interesting, not 
only because it is the spacetime we notice in our day-to-day lives but also because it is 
the marginal case where the phase space of the kinematics is the most constrained, yet 
nondegenerate. We will see the consequences of this below. 

In four dimensions, we use the standard spinor helicity notation (see e.g. [HI 155] ) 
to solve kinematic constraints of the three-point vertex [37]; the momentum bispinor 
form p = XpXp and q = X q X q ensures the null conditions, and the orthogonality condition 
becomes 

= 2p ■ q = (pq) [qp] . (15) 

For real momenta in Minkowski signature, X p and X p are complex conjugates of each 
other (up to a sign). Hence, if (pq) vanishes then [pq] must also vanish and the momenta 
are collinear. If the momenta are taken to be complex, however, the two spinors X p and 
X p are independent. This gives two independent solutions to (1T51) . 

(pq) = or \pq] = 0, (16) 

implying that X p and X q are proportional or, in the latter case, that X p and X q are 
proportional. Since r = A r A r = — X P X P — X q X q , it then follows that there are overall 
two possible solutions: all As proportional, and hence all (ij) vanishing, or all As 
proportional, and hence all [ij] vanishing. We conclude that in D = 4, the momentum- 
space support of the three-vertex exhibits a remarkable twofold branching or chirality. 
We will denote the vertices that are supported by A kinematics by (+) and the vertices 
supported by A kinematics by (— ), as illustrated in figure ED 

The chiral branching implies that the on-shell three-point amplitudes of any 
massless theory in D = 4 have the following kinematic dependence: 

4 +) = 4,(Ai, A 2 , A 3 ) = a (1 2) b (2 3) c (1 3) d , 

4"> = A 3 (X h A 2 , As) = a [1 2] b [2 3] c [1 3] d , (17) 

where a, b, c and d are numbers which may be determined by calculation (b, c, d alone 
may be determined by Lorentz symmetry [56]). For example, the pure-gluon amplitudes, 

% Given the two D-dimensional momenta p and q that satisfies (Tl3|) . we introduce a massive reference 
vector v such that v ■ q = and v -p ^ 0, which can always be found when p and q are not collinear. By 
a Lorentz transformation, we may bring v and p to the forms = (•, 0, . . . , 0) and p^ = (•, •, 0, . . . , 0). 
From the conditions v ■ q — and p ■ q = 0, we conclude that in this frame q^ = (0, 0, •,»,..., •). 
Next we perform a Lorentz transformation in all but the first two spacetime directions, which leaves p 
unaffected and brings q to the form q^ = (0, 0, •, •, 0, . . . , 0). 
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Figure 1. In four dimensions, the kinematics of the on-shell three-point vertex comes 
in two chiral classes: the (+) vertex is supported on As (thus all As are parallel and 
their [ij] vanish), and the (— ) vertex is supported on As (thus all As are parallel and 
their (ij) vanish). 




Figure 2. In a D = 4 cut, one can encounter two three- vertices of the same chirality 
joined by a common line (e.g. see cuts in figure |4jb) and[5])- For the purpose of solving 
the kinematics, one can treat this configuration as an effective four-point vertex of 
definite chirality (note that a corresponding amplitude is not defined). All intermediate 
channels of this effective kinematic vertex are on-shell, as all legs are proportional to 
the same spinor. 



Note that the A 3 cc (l + ,2~,3~) amplitude vanishes for the kinematics of a (— ) 
vertex, as the numerator vanishes faster than the denominator ~ 3 /0 2 . Similarly, 
the A^' tTCC (l~ , 2 + , 3 + ) amplitude vanishes on the (+) vertex. The same is true for all 
three-particle amplitudes in massless gauge and gravity theories, as these amplitudes 
are completely local objects (the spinor-helicity formalism obscures this fact). 

As we are interested in cuts that involve more than one on-shell three vertex, we 
should also understand the properties of configurations of several such vertices. One 
immediate simple "transitivity" property is: if two three-vertices of the same chirality 
are connected by a line with momentum AA, then the two vertices are supported on the 
same chiral branch; thus, all momenta of the two vertices are proportional to a single 
spinor: A or A depending on the chirality. As illustrated in figure El this implies that one 
can treat the kinematic configuration as an effective chiral four-point vertex. Similarly, 
if one encounters a connected chain of three-point vertices of the same chirality, one 
may treat it as an effective chiral higher-multiplicity vertex. All intermediate channels 
of such an effective kinematic vertex are on-shell, as all legs are proportional to the same 



where leg 1 has opposite helicity to that of leg 2 and 3, are 
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(a) (b) 

Figure 3. A quadruple cut of a four-point amplitude in D = 4, and a penta-cut of a 
five-point amplitude in D > 4. In D = 4, we use (+) and (— ) labels on three-point 
vertices to specify the chirality, but in D > 4 the chirality is not a well-defined notion. 




(a) (b) 

Figure 4. A hepta-cut of a four-point amplitude, and an octa-cut of a five-point 
amplitude, both in D = 4. 

spinor (note that there exists no known corresponding effective amplitude, the vertex 
here is used only as a kinematic tool). 

In figures [21 IH and E] we illustrate typical cuts that are built out of only on-shell 
three point vertices. The four-particle quadruple cut in figure [3(a) uses four-dimensional 
momenta and must be built with vertices of alternating chirality [36| 137] . (The non- 
alternating configurations (not illustrated) correspond to degenerate situations where 
the momenta of two external legs are "half-collinear" ; that is, they have a common 
spinor, similar to the situation in figure EJ) In D > 4, the penta-cut in figure [3(b) is 
well defined, but unlike the D = 4 case, the three vertices no longer come with definite 
chirality; therefore, we simply mark the vertices by a light shade (to avoid confusing 
them with tiny loops). In figure HI we show examples of two-loop cuts in D = 4. 
Because the hepta-cut (a) enforces seven on-shell propagators If — 0, one out of the 
eight degrees of freedom of the loop momenta is unconstrained (see discussion below). 
For the D = 4 octa-cut of the five-point amplitude, figure [4(b), the number of on-shell 
conditions matches that of the number of loop momentum parameters; thus, the loop 
momenta is frozen by the cut, similar to the situation of the one-loop quadruple cut. 
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(a) 



■e 



^0 — O — 

(b) 



Figure 5. A non-planar hepta-cut of a two-loop four-point amplitude, and a non- 
planar deca-cut of a three-loop amplitude, both in D = 4. 



Before we conclude this section by working out a concrete example of a M = 4 sYM 
maximal cut, it is convenient to recall some simplifying properties of state sums in 
maximal cuts [3]. If a gluon of helicity (+1) is attached to a chiral (+) vertex, then the 
corresponding fermion and scalar amplitudes vanish, only the pure gluon amplitude has 
support: 



A 



(+),tree/ 1 + 



l 3 V^) 2 /) 3 /; 

4+)^e ( l+ 2" 37)^0. (19) 

The same is true for a gluon of helicity (—1) that is attached to a chiral (— ) vertex. It is 
easy to understand these vanishings from the point of view of charge conservation. For 
a massless fermion the helicity has to be conserved in gluon-fermion interactions; thus, 
in our all out-going notation the helicity must flip. For example, v43 + ^' trec (l+, 27, 37) is 
forbidden by helicity conservation. The case A3 tree (l+, 2^ , 37) is allowed by helicity 
conservation, but it vanishes by virtue of having the wrong chirality. Similarly for a 
complex scalar its U(l) charge has to be conserved in gluon-scalar interactions, implying 







that only amplitudes of the form A 



(+),tree . 



1+ 2+,3 s ) are valid (where in this case ± 



"scalar helicity" refers to the scalar and its complex conjugate), but again this amplitude 
vanishes because of the chirality mismatch. The only nonvanishing amplitude is thus 



A 



(+),tree/-, + 



1+ 2,3), which has the right chirality. This "projection property" allows for 



some great simplifications of the maximal cuts by fine tuning the external states, such 
that all, or almost all, of the complicated particle and state sums in the loops drop 
out. In fact, (TI9~|) is generic to any massless gauge theory (and gravity); however, the 
most useful application is to the M = 4 sYM theory where the amplitudes assemble into 
superamplitudes that are relatively insensitive to the choice of external states. 

If we classify the various state-sum cut contributions according to the helicity 
configuration of the particles (thus ignoring the particle species labels), then the greatest 
simplification that can arise in a nontrivial kinematic solution is the restriction to a single 
allowed helicity configuration for the internal lines. We will refer to this configuration as 
a "singlet". In this case only gluons can propagate inside the diagram, as in figure |6]^a). 
Fermions or scalars are not allowed because in each loop there is at least one vertex 
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Figure 6. A singlet hepta-cut (a) and one of the two helicity configurations (b) in 
a non-singlet cut. The latter allows gluons, fcrmions and scalars to propagate in the 
loop indicated by a dashed circle. The other configuration is obtained by flipping all 
the helicity signs of the lines in this loop. All lines in this figure are cut and carry 
on-shell momenta. The arrows in (a) refer to the direction of momentum flow. 

of the wrong chirality type. The second-strongest simplification allows two helicity 
configurations. In such configurations, the particle content is purely gluonic except for 
one loop in which any particle type can propagate, as shown in figure Mjo). (A single 
fermionic loop always allows two helicity assignments, corresponding to interchanging 
fermion and antifermion, and the same is true for a loop of complex scalars.) 

Now we continue with an explicit example of a maximal cut with only cubic 
vertices. Consider the two-loop four-point hepta-cut in the singlet configuration shown 
in figure [6]^a). Although for analytical computations is not strictly necessary to solve 
the on-shell cut conditions, for illustrative purposes, we do so here. The seven on-shell 
conditions If = 0, together with momentum conservation at the six cubic vertex give a 
set of equations with the solution 

h — £^i^2 > h — h ~ k 2 , <3 = l\ + h\ , 
U = ^4^3^f^ , h — ^3 + U j l& = U — , h = h — h i 
where £ is an arbitrary complex parameter, corresponding to the remaining degree of 
freedom in the integration not frozen by the hepta-cut. The procedure to obtain this 
solution is quite simple. First we note that l\ has to be proportional to the external 
spinors of the vertices that it connects to, namely, l\ oc A1A2. We may try to find the 
overall scale of this momentum, but we observe that there is one degree of freedom 
unfixed by this cut, therefore, we may account for this by introducing an arbitrary 
complex scale £ as in O20p . Thus, we consider l± solved. Next we note that both 
l 2 and I3 are on shell simply by momentum conservation. Therefore, we proceed by 
solving £4. Similar to l\ it has to take the form Z 4 = £'A 4 A3. Again by momentum 
conservation, the neighboring momenta I5 and Iq become on-shell using this ansatz. 
However, the momentum Z 7 is off shell for generic values of £ and Indeed the condition 
= lj = (k^ + U — ^) 2 is the equation that allows us to solve for one of these parameters, 
say The reader may check that this equation is linear in £' and that the solution 
takes the form given in ( 1201 ). 
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The hepta-cut of figure [SJa), for external gluons, is given by 



a (2) _ -7 a tree /(tree /(tree /(tree /(tree /itrei 

A 4 - « ^(1) ^(2) ^(3) ^(4) ^(5) ^(6) 



% 1] [03 <Z X 2> (2Z 2 ) (/ 5 3) (3/ 4 ) 
[/ 6 / 4 ] 3 [/ 2 / 5 ] 3 (Me) 3 



[Ms] 3 (M 2 > 3 (/ 5 / 4 ) 3 



[/ 6 4] [4/ 4 ] [/ 2 / 7 ] [/ 7 / 5 ] (/ 3 / 7 ) (Me) 




3 + ,4~), (21) 



where A^ schematically denotes the three-point subamplitudes of figure [SJ(a), and 
s — (fci + A^) 2 ; t = (&2 + /C3) 2 are the familiar Mandelstam variables. The intermediate 
steps in this calculation can be checked numerically using the solution (|20|) . or by simply 
massaging the analytical expressions using basic properties of spinor products. (For 
example, in order to simplify seemingly complicated products of tree amplitudes, it 
may be useful to keep in mind that maximal cuts are always local objects in the loop 
momenta.) The reader may confirm that (|2~T|) holds for any (nonzero) value of the 
arbitrary parameter £. 

As (1211) is a unitarity cut in the M = 4 theory, we should have summed over all 
particle states as well as helicities in the loops; however, as argued this particular 
kinematic configuration, with conveniently chosen external helicities, is a singlet cut. 
There is only a single purely gluonic contribution to the state sum! For a complete 
treatment of the two- loop four-point hepta-cut, one should also consider other kinematic 
solutions corresponding to other choices of chirality configurations for the cubic vertices. 
In general, this will result in non-singlet cuts, thus requiring the full supersymmetric 
multiplet of states. For such cases, it is useful to employ the four-dimensional on- 
shell superspace formalism of [57J [581 EHJ HSl El [60]; for additional details of on-shell 
superspace in multiloop cuts, see [ST]. However, for the purpose of this review, it is 
sufficient to note that the cuts of the other kinematic solutions give the same result as 
(12T|) . as shown by 

We have now computed the hepta-cut in one of the possible channels of the planar 
two-loop four-point amplitude; this allows us to write up an ansatz for the Feynman-like 
diagram in this "horizontal" channel: a double box scalar integral I^ p \s,t) times the 
coefficient we computed in the cut — s 2 tA i £ ee (l~ , 2 + , 3 + , 4~) . There is a second possible 
channel corresponding a "vertical" double box integral I^ p \t, s), see figured The hepta- 
cut corresponding to this diagram can be computed the same way as the first hepta-cut 
( 12T|) . The result, which we will simply quote, is remarkably similar to ( 12T|) . the only 
difference is a exchange s t (note that the helicity amplitude A\ (1~, 2 + , 3 + , 4~) has 
no cyclic symmetry or flip symmetry that automatically enforces this). Finally, with 
the results of the maximal cuts in our hands we can write down an ansatz for the full 
planar amplitude: 



4 2) (1,2,3,4) = -Af c (l,2,3 



4) {sHl^\s,t) + st 2 I^\t,s 



)) 



(22) 
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Figure 7. An ansatz for the planar four-point two-loop Af = 4 sYM amplitude in 
terms of the scalar integrals I^ F > (s, t) and 1^ (t, s) , as given by the computation of the 
maximal cuts. An overall factor of — A^ ce (l, 2, 3, 4) is suppressed. Other cuts, such as 
near maximal cuts, are needed to detect any terms missing from this ansatz. (However, 
no corrections are to be found as indeed this is the full planar amplitude 



where we have suppressed the helicity assignment (1~, 2 + , 3 + , 4~) that we have been 
considering, the coupling constant, and the usual planar color structure of adjoint 
particles. If, rather, we were to compute the maximal cuts of the remaining four point 
amplitudes, for other gluon helicity assignments (1 , 2 ± , 3^, 4 1 * 1 ), or for the external 
fermions and scalars, the result is the same as (|22|) (this fact can be attributed to 
special properties of the four-point function of M = 4 sYM which hold in all dimensions 
D < 10, see section @J. 

3.4. On- shell three-vertices in higher dimensions. 

In D > 4, the structure of a single three-point vertex is quite similar to the D = 4 
case, as the three momenta p, q, and r can be embedded in a four-dimensional subspace 
(see ( 114]) ) where the details are the same as above. However, the two-fold branching 
is no longer a meaningful notion in D > 4 as one may smoothly connect the two 
chiral solutions of a given subspace by a rotation into the dimensions orthogonal to 
the subspace. Also the simple D = 4 forms of the amplitudes cannot be used in four- 
dimensional subspaces, as the D = 4 spinors (and amplitudes) transform nontrivially 
under the full Z)-dimensional Lorentz rotations. Instead, for the D- dimensional vertices 
one may simply use the fully covariant form of the amplitude, e.g. for three gluons, 

Af»{l, 2, 3) = e x ■ e 2 e^h - k 2 )^ + cyclic , (23) 

which is simply the familiar Feynman gauge vertex contracted with the gluon 
polarizations. Alternatively one may find useful the explicit generalizations of the 
spinor-helicity notation in D = 6 (62} [631 El] an d in D = 10 |64j . 

Since much of the properties of the D = 4 three-point vertex carries over to higher 
dimensions, it should also be true that there are particular choices of gluon polarizations 
that will remove the contributions of the fermion and scalar amplitudes. Indeed, 
choosing the polarization vector of the gluon to be proportional to the momentum 
of one of the neighboring lines is the correct D- dimensional generalization of (I19j) . 

Af e (l 9 , 2 f , 3 f ) = At ee (l g , 2 S , 3.) = , (el oc 
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(a) (b) 

Figure 8. A quadruple cut in D = 4 is a maximal cut, whereas a quadruple of 
the same five-point amplitude in D > 4 is a near-maximal cut, or more precisely a 
next-to-maximal cut. 

Af e (l ff ,2 g ,3 g )^0. (24) 

This is easily seen to be true from the fermion vertex -0 3 e^7 A£ -0 2 , because of the Dirac 
equation k^j^^ = , and from the scalar vertex e^{ki — ^3)/^ because ki ■ kj — for 
three-point kinematics. 

It should be noted that in some cases it is convenient to carry out the method 
of maximal cuts (or other schemes of generalized unitarity) in four dimensions where 
the amplitudes entering the cuts are quite simple and, after arriving at a well-behaved 
ansatz, verify completeness using higher-dimensional spanning cuts. In particular, for 
four-particle amplitudes in the Af = 4 sYM theory, this strategy has been successfully 
carried out through four loops [28j [30] . 



3.5. Cuts with higher-multiplicity vertices 

Maximal cuts do not only include generalized cuts with three-point vertices; according 
to the definition fl9]), higher multiplicity vertices are allowed. This situation occurs, 
for example, in the five-point one-loop quadruple cut in four dimensions, displayed 
in figure Ufa). The four-point vertex may contain a propagator corresponding to a 
pentagon integral, but the constraints of four- dimensional spacetime prohibits all five 
propagators to be on-shell simultaneously (for non-degenerate external momenta). The 
number of cut conditions is saturated by the one-loop quadruple cut in four dimensions, 
and thus quadruple cuts are considered maximal in this spacetime. However, for D > 4, 
the five-point one-loop amplitude evaluated on the penta-cut exists and is considered 
a maximal cut, and the quadruple cut, displayed in figure E^b), would be considered a 
next-to-maximal cut, as per the discussion in section |3~T1 

From a technical point of view, the process of carrying out maximal cuts and 
near-maximal cuts that involve higher-multiplicity subamplitudes is quite similar to 
the case of maximal cuts involving only three-point vertices. The main distinction is 
that higher-point vertices have no definite chirality, not even in four dimensions. This 
implies that the phase space of a four dimensional cut with many higher point vertices 
is in general smoother than four-dimensional cuts with only three-point vertices. The 
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chiral branching of the four dimensional kinematics entering a cut is entirely due to the 
presence of four-dimensional on-shell thee-point vertices. 

For an example of a near-maximal cut with a four-point subamplitude, we will 
consider the two-loop four-point amplitude again. We wish to show that the ansatz 
(122j) receives no more corrections at the level of the next-to-maximal hexa-cuts. First 
we list the different cut topologies that can arise at the next-to-maximal level. Starting 
from the maximal cut topology in figure [61 we may relax the cut condition of any one 
of the on-shell lines. This gives the three cut topologies in shown in figure [9j As before, 
these cuts have to be considered for all possible kinematic solutions (assignment of the 
vertex chirality) and for all possible channels (cyclic and anticyclic permutations of the 
external states). We will work out the figure [9](a) hexa-cut explicitly. Note that the 
vertex chirality and labeling is the same as for figure EJ^a) except that momentum Z7 is 
no longer on-shell. We may recycle the kinematic solution ( 120]) . with the only difference 
that £' is now arbitrary, and thus l-j is off-shell: 

h — £AiA 2 , h — h ~ k 2 , Is = h + ki , 

Z 4 = £,'^4,^3 , / 5 = /c 3 + Z 4 , Z 6 = Z 4 — &4 , l 7 = Z 5 — Z 2 • 

The hexa-cut is given by 



A? 



_ -6 /(tree a tree a tree a tree a tree 
6-cut~ (1) (2) {3) (4) (50 

[M 3 ] 3 (M 2 > 3 (kk) 3 [kkf (hk) 



[h 1] [1 Z 3 ] (Zi 2) (2 Z 2 ) (Z 5 3) (3 Z 4 ) [Z 6 4] [4 Z 4 ] (Z 3 Z 2 > (Z 2 Z 5 ) (Z 5 Z 6 ) 
is 2 t- 



2 1 (14) 4 



'Z? (12) (23) (34) (41) 
= - S 2 4Af c (l",2 + ,3 + ,4-), (26) 

where again A 1 ^ 13 schematically denotes the three- and four-point subamplitudes of the 
cut, in this case figure [9]Ja). The result is surprisingly similar to the hepta-cut ff2T]) . the 
only difference is the appearance of the propagator 1/Z|. This a positive result showing 
that our ansatz for the amplitude passed the test. Indeed the reader may check that 
the same hexa-cut on the scalar double box integral in our convention ([7]) is 

I [P) (s,t) =1 (27) 

6— cut lj 

explaining the appearance of this propagator factor. Finally, we note that evaluating 
all the hexa-cuts of figure |9] for the various kinematic solutions and channels, informs us 
that the ansatz (1221) receives no more corrections at the next-to-maximal cut level. This 
is a strong indication that the ansatz for the amplitude may be complete. To ensure 
that the ansatz is correct, one would have to evaluate it on a complete set of cuts in 
D dimensions, such as the two- and three-particle unitarity cuts (done in ref. [26] in 
D = 4). This we will not do explicitly here (although the next section discusses the 
two-particle cut); we simply note that (1221) is indeed the full planar amplitude [26] . 
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Figure 9. Planar next-to-maximal cuts of the four-point two-loop amplitude. 
4. Some recursive D- dimensional cuts 

In the previous section, we discussed unitarity cuts which rely on cutting all or almost 
all propagators of every internal loop. As a contrast, here we will consider cuts where 
many of the propagators are left uncut. Specifically, will review two types of very useful 
cuts, particular to four-point J\f = 4 sYM, namely the well-known "two-particle cut" 
of the four-point amplitude [26] 133] , and the more recent "box cut" |30j ; both cuts can 
be expressed very simply, and in all generality, in D dimensions, in terms of lower- 
loop expressions. These cuts are useful for computing contributions to the D = 4 — 2e 
dimensionally regulated amplitudes; in this case the regularization scheme implied is the 
four-dimensional-helicity (FDH) version of dimensional reduction [65]. And similarly, 
for amplitudes in D < 10, the corresponding theory implied is the dimensional reduction 
of the ten dimensional M = 1 sYM theory. 

Before introducing the two cuts, we will review a structure present in the four-point 
amplitudes of Af = 4 sYM, which plays an important role in the subsequent discussion, 
due to the especially simple dependence on the external states. It is known that all 
M = 4 sYM four-point amplitudes can be expressed in a common factorized form [30] 

4 L) (1, 2, 3, 4) = g 2L+2 JC(l, 2, 3, 4) [7^(1, 2, 3, 4) , (28) 

where Aa represents the full color-dressed D < 10-dimensional amplitude (c/. the 
color-stripped partial amplitude A4), the kinematic prefactor /C is the same in all loop- 
diagrams, and the loop-specific factor U is independent of the states of the external 
particles. The D- dimensional M = 4 sYM amplitudes are defined as amplitudes of the 
ten-dimensional M = 1 sYM theory dimensionally reduced to D dimensions. All of the 
J\f = 4 state dependence in A4 (16 4 configurations of states) is carried by the kinematic 
prefactor 

/C(l, 2, 3, 4) ee Sl2 S23 AT 2, 3, 4) , (29) 

where v44 rec (l, 2, 3, 4) is the canonically color-ordered D- dimensional tree amplitude. 
Note that multiplying the tree by S12S23 = (&i + k 2 ) 2 (k2 + k^) 2 renders the resulting 
expression fully crossing symmetric; hence, the ordering of the arguments of /C indicates 
no specific color order. 
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In four dimensions, a compact form of K, can be written down using anti-commuting 
parameters r\ in the four-dimensional on-shell superspace formalism [57 1 |58 | 1591 H~9| 15] 160] . 

_ iS^ (Xirji + X 2 r]2 + A 3 ?73 + A 4 7/ 4 ) 

D=4 



A\ ree (l 



2,3,4) 
/C(l, 2,3,4) 



(12) (2 3) (3 4) (41) 
[12] [3 4] 



(30) 



5 (8) (Air/i + A 2 r/ 2 + A 3 r/ 3 + A 4 r/ 4 ) • (31) 

It is not difficult to verify that ([1 2] [3 4])/((l 2) (3 4)) is symmetric under exchange 
of any two legs. Similar compact superspace expressions for the four-point amplitude 
exist in D = 6 [63] using generalized spinor-helicity notation [62], and in D = 10 [66] 
using Berkovits' pure spinor formalism [57]. Related to this, K, represents the color- 
stripped four-point (linearized) matrix elements of the local operator TrF 4 , plus its 
supersymmetric partners. Therefore, /C is a natural prefactor to extract from the four- 
point amplitude in M = 4 sYM. 

All of the color dependence of the amplitude is carried by the external-state- 
independent "universal factor" U^ L \ This factor at tree level, U^°\ is given by 



f/ (0) (l,2,3,4) 



(32) 



s 12 s 23 S12S13 

In general, the universal factor is a sum of L-loop integrals. The integrands entering 
are rational functions of momentum invariants involving the loop and external 
momenta. The universal factors for L = 1,2,3,4, including planar and non-planar 
contributions, may be found by in [3D], and planar L — 5 in [3]. 

The four-point factorization ( )28|) makes it possible to immediately write down 
simple D-dimensional expression for cuts involving this amplitude, usually in terms 
of lower-loop numerator factors. We will now review two related special classes of cuts, 
which owe their simple structure to this property. These cuts can be applied to both 
planar and non-planar amplitudes, making them especially powerful. 



4-1. Two-particle cuts of the four-point M = 4 sYM amplitude 
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Figure 10. A two-particle cut that may be used to construct contributions to 
the {L + L' + l)-loop amplitude from those at L and L' loop orders. 



A regular two-particle cut of a multiloop four-point amplitude has the form shown 
in figure [TD] — it cuts the amplitude into two lower-loop four-point amplitudes. Using 
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the factorization (1281) . it will be possible to write down the cut in terms of lower-loop 
universal factors, thus straightforwardly determining a class of integrals in the amplitude 
that are detectable in this cut. The result of this exercise will be valid in D dimensions, 
whenever the lower-loop universal factors are valid in D dimensions. 

We evaluate the generic two-particle color-dressed cut depicted in figure [101 It cuts 
the (L + L' + l)-loop amplitude A^ +L (1, 2, 3, 4) into the two four-point amplitudes 
■^{ J \— h, 1) 2, l 2 ) and A^ (— Z 2 , 3, 4, Zi) of loop orders L and L', 

4 L+i ' +1) (l,2,3,4) =i 2 V A{ L) (-h, 1, 2, l 2 )Af'\-h, 3,4, /i) , (33) 

2-cut 

AA=4 
states 

where the state sum is over the particles with momenta Zi and Z 2 . 

Using the factorization (1281) and the state- independence of , we can immediately 
rewrite the cut as follows: 

/C(l, 2, 3, 4) f/( i+L ' +1 )(l, 2, 3,4) _ = * 2 C/ (L) Hi, 1, 2, Z 2 ) ^(-'2, 3, 4, Z x ) 

^/C(-Z 1 ,1,2,Z 2 )/CH 2 ,3,4,Z 1 ). 



2-cut 

X 



AA=4 

states 



(34) 



To evaluate this, we use the sewing relation between two four-point color-ordered 
N = 4 sYM trees CUES], 



Y,Af\-h, 1, 2, Z 2 )4 0) (-^2, 3, 4, Z0 =-zsi2S 23 4 0) (1, 2, 3, 4)—^ .(35) 



AA=4 
states 



This sewing relation is valid in any dimension D and for any external states in the M = 4 
multiplet. A straightforward way to confirm (1351) is to work in D = 10 and evaluate the 
sum over states in components, using the fact that in D = 10 M = 4 sYM is equivalent 
to an M = 1 theory composed of a gluon and a gluino. By dimensional reduction the 
sewing relation fl35l) then holds in any dimension D < 10. Recently, this equation has 
also been verified directly in six dimensions using an on-shell superspace [53] . 

Substituting in the relation A± = /C/(si2S23) given in (|28l) . we find that (135]) is 
equivalent to 

J2 £Hi, 1, 2, Z 2 ) /C(-Z 2 , 3, 4, Zx) = - i s\ 2 JC(1, 2, 3, 4) . (36) 



states 



Applying ff36j) to (]34p . we find the key equation for building all contributions from 
two-particle cuts directly in terms of the Us: 

U {L+L ' +1 \l, 2, 3, 4) = i s\ 2 U^(-h, 1, 2, Z 2 ) U^ L '\-l 2 , 3, 4, Z x ) . (37) 

2-cut 

Equation (1371) is rather powerful. No complicated calculations remain in order to obtain 
all contributions visible in two-particle cuts; they are given simply by taking the product 
of lower-loop results. The color-dressed £/( L+L ' +1 ) is given immediately as a sum over 
products of individual integrals residing inside the and U^ L '^ factors, up to terms 
that vanish because of the on-shell conditions, — = 0. 
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Figure 11. Sample contributions to the full color-dressed two-particle cut for 
L = 2 and L' = 1. The diagrams on the right-hand side show some of the 
terms in f/( 4 )(l,2,3,4) that are constructed from these cuts, using and 
U^ l \ The explicit color factors, as well as factors of i, have been omitted. 



P 

> i(p + q) 2 



<1 

Figure 12. The rung rule |26[ [33] for generating a higher- loop contribution by 
inserting vertical rung between two lines inside a loop diagram of N = 4 sYM 
amplitudes. The dimension of the extra propagator is compensated by the 
extra factor i(p + q) 2 in the numerator. This rule reproduces the numerators 
of planar diagrams generated by iterated two-particle cuts. 



Figure [TT] illustrates diagrammatically some of the terms generated by ( )57j) for 
the case L = 2 and V = 1. For simplicity, we draw only the planar contributions of 
U^'(—li, 1, 2, l 2 ), omitting the color factor f abc s as well as all factors of i (which can be 
restored at one's convenience). The denominator factors in f/( 2 ) and correspond to 
the un-cut propagators that are graphically visible on the left- and right-hand sides of 
the dashed line, respectively. Similarly, the numerator factor for each parent graph on 
the right-hand side of figure [11] is given by forming the product of the numerator factors 
for the one-loop box diagram and two-loop double box (taking into account the proper 
permutation of legs), and then multiplying by two powers of s\2- 

The simple formula ( 157)) for the two-particle cuts provides a rather useful tool 
for generating many higher-loop contributions by recycling lower-loop ones. Indeed, 
it is the mechanism behind the many of the diagrams given by the well-known "rung 
rule" [26] [55] , shown in figure [12] Unlike the two-particle cut formula, the rung rule is 
a heuristic rule that should be applied with some care; it can potentially suggest non- 
contributing terms when the two-particle cut does not apply [3j. That warning aside, for 
planar M = 4 sYM, the rung rule is rather powerful. For example, it is quite interesting 
that the rung rule generates terms that are invariant under the dual conformal symmetry 
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Figure 13. The box-substitution rule [I] for generating a higher-loop 
contribution by inserting a one-loop four-point box subintegral into a four- 
point vertex. In this example, we substitute a box into the central four-point 
vertex in the four-loop "window" diagram. The result is a five-loop integral 
that cannot be obtained from two-particle cuts or the rung rule. 



of M = 4 sYM [3j H], presaging modern awareness of this symmetry. 
4-2. Box cuts of multiparticle M = 4 sYM amplitudes 

The two-particle cut has a particularly nice generalization that also relies on the simple 
structure of the four-point amplitude of M = 4 sYM. This cut can be thought of as 
a (generalized) four-particle cut that isolates a four-point subamplitude, which in the 
simplest nontrivial case is a box subdiagram; thus, the cut was named "box cut" in [30J. 
A crude version of this cut appeared already [1] as a "box-substitution rule" . It allowed 
the construction of L-loop contributions with a box subgraph, starting from (L — l)-loop 
contributions with a contact interaction, as illustrated in figure [13j Similar rules were 
discussed in conjunction with the leading singularity method in [53]. Unlike the two- 
particle cut discussed in the previous section, the box cut is applicable to multiparticle 
amplitudes. Now we will review the derivation of the L>-dimensional box cut. 
Consider the generalized cut of an L-loop n-point amplitude, 

A {L) 



box cut ; 

Af=4 



that is composed of a generic set of color-dressed amplitude factors, except for the i th 

.(£') 



such factor, which we take to be a color-dressed L'-loop four-point subamplitude, A, 



Example of such box cuts are given in figure [TH 

As mentioned, the L'-loop four-point subamplitude of M = 4 sYM is special because 
of the factorization property fl28|) . Labeling the cut momenta by l\, I2, 13, U, we have 

A[ L ^(h, Z 2 , l 3 , U) = 4% (h, Z 2 , l 3 , k) (h+kfik+h) 2 U<- L \h, h, h, k) , (39) 

where, as in the previous section, we use A4 to represent the color-dressed amplitude, 
and only the color-ordered tree amplitude factor A^L depends on the states crossing 
the cuts. Therefore, we can pull the factor out of the sum over states in (138j) . 

leading to a simpler expression in the summand: 

A L) j box = {k+hfih+hfU^ (k, l 2 , i S) k) A (D ■ ■ ■ • • ■ Am) ■ (40) 

cut J\f=A 
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Figure 14. Two examples of multiparticle multiloop "box cuts" . They reduce 
to lower-loop cuts by replacing the red four-point subamplitudes by four-point 
color-ordered trees, multiplied by known numerator and denominator factors. 
This property allows these cuts to be computed easily in D dimensions, once the 
amplitudes with fewer loops are known. White holes represent loops and the 
darker (red) subamplitudes mark four-point amplitudes amenable to reduction. 



The state-sum is identical to a lower-loop cut, that of the (L — L')-loop amplitude, but 
utilizing the color-ordered contribution to the i th tree. This fact immediately gives a 
simple relation between the L-loop box cut and contributions to the reduced (L — L')- 
loop cut under the same cut conditions. 

We can formally write down an equation relating the cut of an L-loop amplitude 
to a cut of a lower-loop one as, 

A^ = (h + h) 2 (l2 + h) 2 U( L '\h,l 2 ,l 3 ,k)A^ . (41) 

box cut cut 

We introduced the reduced cut A notation to emphasize that the state-sum in (1401) 
is exactly a (L — U) loop unitarity cut which is color-dressed with f abc everywhere, 
except for the four-point color-ordered tree amplitude whose associated color factors 
are accounted for in the L'-loop universal factor U^ L '\ 

Figure [15] gives an example illustrating that the box cut is very a practical 
tool to work with. Using a back-of-the-envelope calculation, one can determine the 
numerator polynomial of a four- loop integral [31 [30] , using the known three- loop four- 
point amplitude [261 I2H]- Although this particular example is planar, it is just as 
simple to use the box cut for non-planar contributions. For example, inserting the box 
into the four-point vertex in figure [15] in a non-planar fashion generates non-planar 
four-loop integrals. 

The box cut is an extremely efficient way to obtain contributions to parent graphs 
that contain a lower-loop four-point subgraph. As mentioned earlier, the LI — 1 box cut 
is closely related to the box-substitution rule depicted in figure [£31 The box cut also 
generates contributions that are consistent with the two-particle cut and/or the rung 
rule [26] . 

From the above covariant derivation, it follows that box cuts are valid in any 
dimension, if both the reduced cut A^t L ^| cut and the four-point universal factor U^ L '^ 
are known in D dimensions. As a practical matter, the universal factors entering the 
lower-loop amplitudes should already be known in D dimensions, prior to attempting 
the higher-loop calculation in D dimensions. 
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Figure 15. Practical application of the box cut to determine the numerator 
N for a four-loop integral. The (red) dashed cut conditions around the upper 
left box in this diagram allow us to replace it by the product of a reduced cut 
diagram, some kinematic factors and a box integral. The reduced cut diagram 
is then expanded into two three-loop "tennis-court" diagrams, corresponding 
to the two allowed channels of the marked four-point vertex. The relevant 
kinematic pieces of the tennis-court diagrams, i.e. the numerators and the 
spurious propagators, are extracted from the known three-loop contribution. 
Assembling all the kinematic factors gives the result for N (with the overall 
factor fC suppressed), which is free of spurious propagators. 



5. Duality between color and kinematic factors 

In modern S-matrix schemes, one successfully constructs multiparticle and multiloop 
amplitudes using minimal input from the theory as defined by the Lagrangian. Rather 
than employing the nuts and bolts of the action, one relies heavily on quantum field 
theoretic consistency. For example, a wealth of tree-level amplitudes for interesting 
theories can be fully constructed through BCFW [HE] using only the on-shell three- 
point amplitude as input. Similarly, as discussed in section [3] and in greater detail 
in [38], the unitarity method allows the construction of all loop level contributions 
using as input only lower-loop subamplitudes. A natural question is whether one 
can further minimize the theory-specific information needed to construct a full color- 
dressed scattering amplitude. The conjectured duality between color and kinematic 
factors [321 12H] described in this section suggests such an approach. 

The duality relies upon the cubic-graph representation of gauge theory amplitudes, 
introduced earlier (j2J), which we write again as 

ML) _ ;L U-2+2L f d LD i 1 NiQ 
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This time we have made the dependence on the coupling constant explicit, and pulled 
out an overall phase relative to fl2]). We note that in general the N cannot be uniquely 
specified; if calculated by Feynman diagrams, they are gauge choice dependent. If 
calculated by other means, they are allowed to transform under a generalized gauge 
transformation [32} [29], EQ], which corresponds to the ability to move contact terms 
between different graphs, while leaving the amplitude invariant. This freedom allows 
for many different representations of the same amplitude, some of which make special 
properties of the theory manifest. 

One very special property is the duality between color and kinematics. This duality 
asserts the existence of representations (142jl where the kinematic numerators Ni obey 
the same general algebraic properties as the color factors Cj. Specifically, the N obey 
Jacobi relations and have antisymmetric cubic vertices analogous to the color factors, 
schematically: 

Ni + Nj + N k = <(=> Ci + Cj + C k = 0, (Jacobi identity) (43) 

N — > — N Ci — > —Ci (vertex-flip antisymmetry) 

where the first line signifies the Jacobi identity valid for specific triplets of graphs in the 
amplitude, and the second line represents the action of flipping the ordering of a cubic 
vertex in a graph. In addition, it is often convenient at loop-level to impose the self- 
symmetries or graph automorphisms of the Ni, similar to the self-symmetries obeyed by 
the color factors Cj. 

The duality between color and kinematics was first proposed by Bern and the 
present authors at tree-level [32], and later conjectured to hold to all loop orders [29] in 
generic gauge theories. It generalizes a four-point property of tree-amplitudes observed 
in cuts of multi- loop-amplitudes [32] , and also noticed in studies of zeros of scattering 
amplitudes in the 1980s |69j . 

Many of the the interesting and surprising consequences of this duality will not be 
treated in detail, as they are outside the scope of this review. Nevertheless, in section l5\4l 
we will briefly comment on some of these consequences, as related to tree amplitudes, 
Lagrangian manifestation, and constructing gravity amplitudes. For the main body of 
this discussion, we will instead focus on how the duality can be used in the explicit 
construction of multiloop gauge theory amplitudes. We will see that the constraints 
imposed by the duality allow us to drastically reduce the necessary input from unitarity 
cuts. (As long as the duality remains a conjecture, it is of course necessary to verify any 
obtained amplitude on a set of spanning cuts.) We may loosely quantify the amount of 
theory-specific information required for the construction of an amplitude by counting the 
number of necessary maximal and near-maximal cuts. Alternatively, one may count the 
graphs used to represent the full amplitude and compare this number to the very small 
number of graph numerators not completely fixed by the kinematic Jacobi relations (I43j) . 
We call the latter ones "master graphs"; from these, all other contributing graph 
numerators can be expressed using linear relations. In section 15.31 the four-point three- 
loop N = 4 sYM amplitude, comprised of 12 graphs, is constructed using a single master 
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Figure 16. A Jacobi identity between three cubic four-point two-loop topologies. 
These three graphs have the same connectivity except for the shaded (pink) edge. 
Every Jacobi identity and kinematic identity between numerator factors can be 
graphically seen as a relation local to a set of three edges. For the Af = 
4 sYM amplitude, only graphs (P) and (NP) contribute [26] . 

graph. The numerator of this single master graph, after imposing all symmetries and 
duality relations, is completely fixed by a single maximal cut, illustrating the dramatic 
simplification involved. 

Nontrivial evidence supporting the duality conjecture at loop level is found by 
studying the four-point amplitude of the Af = 4 sYM theory at one [70], two [26] . 
three [29], and four loops [71], and for the five-point amplitude through three loops [72] . 
Evidence that the duality holds even in nonsupersymmetric theories is provided [29] by 
the two-loop four-gluon identical- helicity pure Yang-Mills amplitude [73] . 

We will elaborate on the duality using graph operators, but first it will be instructive 
to look at an venerable example from the literature where the duality ( 1431) is known to 
hold. 



5.1. Two-loop example 

An instructive nontrivial example can be found using the four-point two-loop M = 4 
sYM amplitude considered previously in section [3j As we have already shown how 
to construct the planar part using other methods, here we merely verify that the full 
amplitude, given in the literature [26], satisfies the color-kinematics duality. In the 
subsequent discussion we will show how such amplitudes can be directly constructed 
using the duality. The two-loop amplitude is 

4 2) = -^ (P) ^34/^ +perms), (44) 

where the two scalar integrals contributing are drawn in figure [161 an d the "perms" 
instructs us to sum over all distinct permutations of external labels. For convenience 
we have pulled the kinematic (and color) factors outside the integral, as they do not 
depend on the loop momenta. The numerators and iV"( NP ) are given by 

N (P) = N (m =JC ( kl + k 2 f , (45) 
where K is the state-dependent prefactor introduced in ( 1291) . 
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The color factors and C^ NP \ obtained by ([3]), are related through the Jacobi 
relation shown in figure HH to a penta-triangle graph C( FT \ 

C (P) = C (NP) + C (PT) (46) 

If ( 144)) satisfies the color-kinematics duality, we should expect a similar relation for the 
kinematic factors: 

iV (P) =iV (NP) +iv (PT)_ (47) 

This indeed holds. The amplitude (jUJ) does not involve any penta-triangle integral; 
therefore, we set A r ^ PT ^ = 0, and thus ( 1471) is consistent with the result in ( )45|) . 



5.£. Graph operators and the color-kinematics duality 

In order to phrase the duality precisely, we will now introduce graph operators. This 
is not an exercise in pedantry. Indeed, precisely defining the duality in terms of graph 
operators allows for natural automation — especially useful when considering higher-loop 
contributions that may involve many graph topologies. 

The color Jacobi identity directly relates graphs that share entirely the same 
connectivity except for one edge, e.g. the shaded (pink) edge in figure [16j The 
language we will use to specify this relation will involve operators that take a graph 
and a particular edge of that graph, and rearrange that edge's connectivity so as to get 
another graph. We can specify any internal edge e of a cubic graph in terms of the four 
edges that touch its bounding nodes. If the order of edges around one node is (a, b, e), 
and the order of edges around the other node is (— e, c, d), we can specify the edge e as 
e = {(a, b), (c, d)}. With this terminology we are ready to consider some mapping t to 
take a graph Q and one of its edges e, to another graph Q t {e)'- 

i(G, e) = g t[e) , (48) 

where Q t (e) has the same connectivity of Q except the nodes (a, b, e) and (— e, c, d) have 
been replaced with (d,a,t) and (—t,b,c), thereby defining the edge t = {(d,a), (b, c)}. 
We will also use 

u(G, e) = £? u(e ) , s(G, e) = Q , (49) 

where Gu(e) has the same connectivity of Q except nodes (a, b, e) and (— e, c, d) have 
been replaced with (d, b, u) and (—it, c, a), with edge u = {(c, a), (d, b)}. The s operator 
simply returns the original graph unchanged independent of the edge specified. These 
operations are depicted visually in figure [13 

We are now ready to concisely specify the color Jacobi identity associated with 
every edge e of every graph Q, 

= c [s(g, e)] - C [t(Q, e)] - C [u(Q, e)] 

= c[g}-c[g t{e) }-c[g u{e) ], (50) 
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Figure 17. Three graph operators taking (graph Q , edge e) to graphs Gt c ,Gu c and 
Q. All three graphs have the same connectivity, except for the edge touching edges 
(a, b, c, d). 



where C[Q] gives the color factor arrived at by dressing every node of graph Q with the 
jabc structure constants. The color-kinematics duality is satisfied in a representation 
where for all edges e of all graphs Q, 

= N [s(g, e)] - N [t(g, e)] - N [u(Q, e)] , (51) 

where N[Q] gives the kinematic numerator factor associated with the graph Q. 

Additionally, we may use this formalism to write down some relations that originate 
from the antisymmetry of the cubic vertices 

C[Q{( d ,a),(b,c)}] = ~ ' 

N[g { ( d>a ),(b,c)}] = - N[g {Mi ( bjC ) } ] . (52) 

In the following section, we will see that the set of constraints on the numerator factors 
f l5Tj) and (|52|) can be understood as a set of functional equations that must hold for the 
duality to be satisfied. 
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Figure 18. Three-loop four-point cubic graphs considered in the main text. The 
external momenta is outgoing and the shaded (pink) edges mark the application of 
kinematic Jacobi relations used in (|54")l . Note that only graphs (a)-(l) contribute to 
the AT = 4 sYM amplitude where the duality between color and kinematics is made 
manifest. 



5.3. Three-loop example 

In this section, we reexamine the four-point three-loop M = 4 sYM amplitude using 
the duality between color and kinematics [29]. This amplitude was originally given in 
[271 EH] in terms of nine cubic diagrams. For this exercise, we start by considering a 
larger set of 25 graphs, which are related to any of the original nine diagrams by a 
single application of a kinematic Jacobi relation. However, eleven of these diagrams 
contain triangle subgraphs, which the no-triangle property of M = 4 sYM [1] suggests 
will not contribute. After removing those with one-loop triangle subgraphs we have the 
14 graphs depicted in figure [181 We will see that this set of diagrams is sufficiently large 
to admit a manifest representation of the duality. 

Now we will introduce the kinematic Jacobi relations that the numerators of each 
diagram must satisfy. Each numerator depends on three independent external momenta 
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and (at most) three independent loop momenta. While in general a graph's numerator 
factor will depend on the states of all external particles, the four-point amplitude of 
M = 4 sYM is special in that all of its state dependence can be captured in the overall 
factor /C described in section HI Since this factor is common to all terms and the 
kinematic Jacobi relations are homogeneous, this factor plays no role in the equations. 
With this in mind, we suppress all state dependence but make the dependence on 
momenta explicit: 

= N^ x \k 1 ,k 2 ,k 3 ,h,k,l 7 ), (53) 

where the canonical labeling for each graph x G {a, b, c, d, e, f , g, h, i, j, k, 1, m, n} is given 
in figure [TBI 

For every internal edge and for every graph there is one functional equation relating 
the numerator of that graph to the numerator of two other graphs through the kinematic 
Jacobi relations. The 14 graphs of figure [18] each have ten propagators, giving 140 
relations in total (before accounting for possible overcounts). One can write down 
many more equations using the graph automorphisms of each numerator - ensuring that 
numerators are invariant under equivalent ways of labeling symmetric graphs. However, 
this large set of equations contains much redundancy; the number of independent 
functional relations is quite small. We will not elaborate on how to methodically 
reduce this system; instead, we observe that the following 13 kinematic Jacobi relations, 
together with all automorphism relations, solve all constraints imposed by the color- 



kinematics duality: 








= N^(k 1 ,k 2 ,k 3 ,l 5 ,k,l 7 ), 
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N W 






AT(c) 




(Jo) 


AT(d) 


= N^(k 3 , h, k 2 , l 7 , l 6 , k 1)3 -l 5 + l 6 - h) 






+ (k 3 , k 2 , k u l 7 , k, fc 2 , 3 + h - h) , 
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N® 


= N^(k 1 ,k 2 ,k 3 ,l 5 ,l 6 ,l 7 ), 


(Jf) 


AT(g) 


— N^^kx, k 2 , k 3 , l 5 , l 6 , h) , 


(J,) 


AT(h) 


= -N<*> {h, k 2 , k 3 , h, / 6 , k 1>2 -h- l 7 ) 
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(Jh) 




= iV( c )(A; 1 ,A ;2 ,A ; 3,/5,/ 7 ,/6) 






- iV (c) (k 3 , k 2 , fci, -k 4 l 6 , -l 6 - h, h) , 


(Ji) 


JVG) 


= N®(k u k 2 , k 3 , l 5 , l 6 , l 7 ) - N^(k 2 , h, k 3 , l 5 , l 6 , l 7 ) , 




AT(k) 


= 7V« (h, k 2 , k 3 , l 5 , k, l 7 ) - N® {k 2 , h, k 3 , l 5 , l 6 , l 7 ) , 
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= 0, 
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where foj = —(k±- 


- k 2 + k 3 ), kij = ki + kj, and for notational simplicity 



(54) 



suppressed the arguments of the numerators on the left-hand-sides; these all have the 
canonical argument (k±, k 2 , k 3 , Z 5 , Z 6 , l 7 ). The 13 relations are labeled by (J x ) which 
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correspond to Jacobi relations specified by the marked edges in figure O In the above 
relations we have immediately set to zero any diagrams that do not feature in figure [18j 
This explains why some relations in ( )54"|) involve fewer than three terms. 

Analyzing the result of (154]) . we immediately see that the color-kinematics duality 
has constrained and to vanish, given that we insist on the absence of graphs 
with one- loop triangles. We see that N^ g \ are directly expressed as linear 

combinations of the function N^ c \ Furthermore, N®, A"W are expressed as linear 
combinations of the functions N^ g \ and thus they can also be solved in 

terms of N^ c \ Finally we note that = = = and that can 

be expressed as a linear combination of . So each of these numerators can also be 
solved in terms of . Therefore, every numerator is entirely expressible in terms 
of the planar diagram numerator 

Although the solution to (1541) is quite constraining, it does not enforce all of the 140 
Jacobi-like relations we began with, nor does it enforce all automorphism relations. As 
mentioned, we will obtain the full solution by enforcing automorphic invariance of each 
graph numerator. Using the partial solution of ( 154"|) . these relations can be translated 
to self-constraints on N^ e \ For example, the figure H8T e) diagram is invariant under 
reflections up -H- down; thus, must satisfy 

(k u k 2 , k 3 , l 5 , l e , l 7 ) = iV (e) (k 2 , h, k 4 , k lf2 - h, l 7 , h) . (55) 

Similarly, figure [TBlj ) is invariant under 1 -H- 2. However, just like color factors, 
numerators are odd under the reordering of a cubic vertex; therefore, we have 

N®(ki, k 2 , k 3 , l 5 , l 6 , h) = -N®(fa, h, k 3 , l 5 , l 6 , l 7 ) . (56) 

In fact, the reader may check that this equation is already enforced by ( 1541) ; thus, 
this particular automorphism relation does not give rise to any new constraints on 
. Beyond these two examples, every diagram but figure [T8T g) is constrained by 
automorphism relations, which through ( 15~4"1) translates to possible further constraints 
on iV( e ). 

One way to efficiently solve the remaining kinematic Jacobi and automorphism 
relations is to write down an ansatz for N^. As mentioned, we assume that the state- 
dependent prefactor is K, = si2S23^4 ree (l, 2, 3, 4). We write the remaining factors as a 
local degree- two polynomial in 17 momentum invariants: 

^inLtz = £ P i S l2, S23, ki ■ lj, li ■ Ij) , (57) 

where ki = ki,k 2 , k 3 and = l 5 ,l 6 ,l 7 , and the degree of the polynomial is given by simple 
dimensional counting. Although we could use various facts about the power counting 
of = 4 sYM to simplify the ansatz, we proceed with the most general homogeneous 
degree-two polynomial of 17 variables, which gives 17 x (17+ l)/2 = 153 parameters to 
be determined. 

Using such an ansatz, we find that only the following four automorphism relations 
are required in order for all automorphic relations to be satisfied: 

iV (c) (h, k 2 , k 3 , l 5 , l 6 , l 7 ) = N {c) (h, k 2 , k 4 , l 5 , l 5 - l 6 , k 3A - h) , 
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iV (f) (k u fc 2 , h, l 5 , Z 6 , l 7 ) = iV (f) (fc 2 , kx, fc 4 , fc 1>2 - Z 5 , h-l 5 - l r , fc 4 + Z 5 - Z 6 ) , 

(fc 1; fc 2 , fcg, / 5 , / 6 , / 7 ) = iV (h) (fcl, fc 4 , fc 3 , fcl - h, l 7 , k) , 

(fci, fc 2 , fc 3 , Z 5 , Z 6 , Z 7 ) = ^ (h) (fc 2 , fcj, fc 4 , fc 1)2 - Z 5 - Z 7 , -Z 6 , Z 7 ) . (58; 

Enforcing (|54|) and (!58j) . the ansatz is constrained from 153 parameters to only three 
free parameters, which we label a,/3,7, giving 



N 



(e) 



sol 



K 



a S12S23 + P (s 2 i2 + S23) + 7 (s? 2 + 2(s i2 - s 23 ) fci • h 
+ 2(2si2 + s 23 ) k 2 -l 5 + 6S12 fc 3 • Z 5 



(59) 



After checking that this also solves the original 140 Jacobi relations and all 
automorphisms, we conclude that (159]) is a threefold family of potential solutions that 
makes the duality between color and kinematics manifest. 

We still need to ensure that these solutions actually describe the N = 4 theory. 
One such source of theory-specific data can be a generalized unitarity cut. We consider 
the maximal cut of graph (e), using the known expression in the literature [26] . 



N (e) 



2/C sio ki • Zk . 



cut 



Then, we compare to the maximal cut of the solution (1591) . 



(e) 



sol 



cut 



a S12S23 + /3 {s\ 2 + S23) + 7 (-512^23 - 6S12 fc 4 ■ Z 5 ) 



-S12 



(60) 

(61) 
2fc 4 • Z5, 



where we have used the relations fci ■ Z5 = 0, 2fc 2 ■ Z 5 = si 2 , and 2fc 3 ■ Z 5 
valid on the maximal cut. 

Comparing (|59|) and (|60|) . we find that a = 1/3, j3 — 0, 7 = —1/3, giving the fully 
solved numerator 



iV (c) = i/C ( Sl2 (s 23 - s 12 ) - 2(s 12 - s 23 ) h ■ Z 5 
- 2(2si2 + s 23 ) fc 2 • k ~ 6s i2 fc 3 • h 



(62) 



All 12 nonvanishing numerators can now be obtained from ( 16 2 p . using (I54p . They are 
explicitly given in [29]. Together, they represent a duality-satisfying expression of the 
three-point four-loop M = 4 sYM amplitude. This was verified on a spanning set of 
cuts in |29|. 



5.4- Tree level, Lagrangians, and gravity 

While this review focuses on multiloop methods for gauge theory amplitudes, we 
should mention a few other interesting results intimately related to the kinematics- 
color duality. Consider the formula ( 1421) in the tree-level limit L = 0. This gives the 
cubic representation of tree amplitudes introduced in [32] : 

jg 73 flxfl-zflz fl m 
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where the sum runs over all cubic tree-level graphs 73 (including all distinct permutations 
of external legs), which all have unit symmetry factors. For notational uniformity 
with [32] we use = iV* ree and q = C* ree . 

As shown in [32], assuming the duality between tree- level color and kinematics 
leads to novel relations between color ordered tree amplitudes. These decompose any 
color-ordered partial amplitude into a sum of (n — 3)! basis amplitudes, schematically, 

ArKl),a(2),...,a(n))= £ K^A^ (1, 2, 3, p(4), . . . , p{n)) , (64) 

peSVi-3 

where a and p are permutations, and Kp are some universal (state-independent) non- 
local kinematic coefficients. We will not give these factors here; they can be deduced 
using the results in [32] . Similar relations between partial amplitudes were subsequently 
found in string theory [74] using monodromy of integration on the worldsheet. In the 
field-theory limit, these string relations provided a first proof of ( J64l) . Using BCFW 
recursion, the same relations were later derived from field theory [75] . In [76] , analogous 
relations was shown to exist in non-commutative Yang-Mills theory. 

As for the duality itself, phrased in terms of color and kinematic numerator 
factors and corresponding Jacobi relations (1431) . less is known. An investigation in 
this direction is [77] where the framework of heterotic string theory is used to shed light 
on the duality. In addition, there are known explicit solutions for tree-level duality- 
satisfying numerators to all multiplicity [7H1 [79] . While these solutions demonstrate 
the existence of numerators satisfying the duality, they take a particularly non-local 
and asymmetric form (under relabeling of external states). The existence of local or 
symmetric forms of duality-satisfying tree-level numerators to all mulitplicity remains 
an open question. 

In [80], it was shown that the duality between color and kinematics could be made 
manifest at the Lagrangian level of Yang-Mills theory, by introducing higher degree 
interactions at the five-point and six-point levels, schematically: 

£ym = -hF^F^ u + £ 5 + £ 6 + . . . , (65) 

and at the same time introducing auxiliary fields so to make all interactions cubic. 
Remarkably the corrections C 5 +£ 6 -\-. . . vanish by the usual gauge-group Jacobi identity, 
and do not alter the theory in any way other than modifying the relative values of 
individual Feynman diagrams. 

Finally, we remark that one can easily construct gravity amplitudes using the 
duality; gravity numerators are double copies of gauge theory kinematic factors |32l [29] . 
Given an n-point L-loop gauge theory amplitude with duality satisfying numerators iVj, 
a gravity amplitude takes the form 

where iV* are a set of kinematic numerators for the amplitude of a possibly different 
gauge theory (which need not explicitly satisfy the duality), and k is the gravitational 
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coupling constant. This formula is the generalization of the analogous tree- level 
(L = 0) formula [32]. Assuming the color-kinematics duality, the double copy property 
was proven [80] recursively in n for the explicit case of maximal supersymmetry 
[A/ — 4] ® [A" = 4] — > [A/" = 8] and similarly for pure Yang-Mills theory going to pure 
gravity plus axion and dilaton. At the classical level, ( )66l) encodes (and arguably 
explains) the field-theory forms of the string-theoretic Kawai-Lewellen-Tye (KLT) tree- 
level relations [81] between gravity scattering amplitudes and gauge theory amplitudes. 
Intriguingly, the double copy form of gravity ( 1661) appears to generalize seamlessly to 
the quantum level, to any loop order [29] . 

We close this section by referring the reader to the literature for further application 
of the duality, the new tree-level amplitude relations, and the gravity double copy 
property jH2J [30j [79j [66] ; see also reviews [83] . 

6. Conclusions 

In this review we have covered some recent methods for constructing compact integral 
representations of multiloop scattering amplitudes, applicable to planar and non-planar 
gauge theory in general, and to A/* = 4 sYM theory in particular, valid in four dimensions 
and in D > 4. 

The maximal cut method is a convenient scheme in the generalized unitary 
framework that allows the construction of amplitudes in any massless theory starting 
from the cuts with the most on-shell propagators. The method has a natural hierarchy 
which determines the order in which unitary cuts should be evaluated in order to capture 
missing contributions of an initial amplitude ansatz. Maximal cuts are followed by 
next-to-maximal cuts and thereafter next-to-next-to-maximal cuts and so on, until the 
full amplitude has been determined. These maximal and near-maximal cuts can be 
evaluated in four dimensions and in D > 4, for planar as well as non-planar cuts. A 
close relative to the maximal cut method is the leading singularity method [531 E] , 
which focuses on the A* = 4 sYM theory in four dimensions. 

Complementary to the maximal cut method, we have reviewed two particularly 
simple .D-dimensional multiloop unitarity cuts in the M = 4 sYM theory: the ordinary 
two-particle cut of the four-point amplitude, and the box cut. These two types of cuts 
have a recursive structure that can be implemented using simple pictorial calculations. 

Although the M = 4 sYM theory naturally lives in four dimensions, the D- 
dimensional cuts and amplitudes are important for several reasons. Since M = 4 sYM 
amplitudes exhibit infrared divergences at loop level they need to be regulated, either 
using dimensional regularization, or using mass regulators [21J. The former method 
naturally requires D-dimensional amplitudes, and the latter can in many cases be 
rephrased in terms of D- dimensional amplitudes, where the higher-dimensional momenta 
effectively act as masses from a four-dimensional perspective. Beyond this, there is, of 
course, the inherent theoretical interest in studying higher- dimensional theories. 

Using the duality between color and kinematic factors we considered how to 
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compute multiloop integrands in a very efficient manner. We discussed the three-loop 
four-point amplitude of M = 4 sYM in some detail, showing how the duality can be 
made manifest using 12 nonvanishing cubic three-loop graphs, all determined from only 
one master graph. The duality, which is conjectured to hold for any multiparticle and 
multiloop amplitude in M = 4 sYM, as well as less-supersymmetric theories (including 
pure Yang-Mills), hints of some deep field-theoretic origin, which once understood 
will explain the suggestive new ways to organize and calculate scattering amplitudes. 
Similarly, the connection to gravity amplitudes as double copies of gauge theories, similar 
to the string-theoretic KLT relations, hints of an underlying unifying framework that 
includes a particularly novel notion of graviton compositeness. From this perspective a 
connection to string theory is not surprising. 
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